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I. INTRODUCTION 

Piecewise-flat (PF) simplicial manifolds are a crucial computational framework for systems 
with dynamic geometry or systems with complex topologies /geometries. In general relativity, PF 
manifolds are a cornerstone of the coordinate-free discretization introduced by Regge [1] often 
referred to as Regge calculus (RC). RC is often regarded as the backbone of the semi-classical or 
low-energy limit of, or even an effective theory [2] of, quantum gravity. It has also proven to be a 
useful computational tool in numerical relativity [3] . For numerical solutions to partial differential 
equations, PF manifolds are one clear method to discretize complex geometries over which the 
differential equations act [1] or in conformal transformations on 2 dimensional surfaces of arbitrary 
genus [5]. 

One may view RC as an approach to characterizing the intrinsic geometry of PF manifolds. 
In other regards the exterior calculus prescribes the calculus of fields on the curved background 
of smooth manifolds. Exterior calculus approaches to functions and fields on and the curvature 
tensors of the PF manifold have been of great use in preserving geometric notions on the discrete 
manifolds. In the canonical approach to RC, one can describe an exact action principle [6, 7], yet 
the discrete curvature tensors give only approximate expressions after smoothing over the disconti- 
nuities. The exterior calculus methods were also quite useful in developing a clearer understanding 
of the Einstein tensor in RC [8, 9]. In the numerical analysis for PDE's on PF manifolds, discrete 
exterior calculus (DEC) is used to preserve geometric symmetries of continuous systems on discrete 
manifolds [10] and for providing the bedrock on which to construct geometric fiows on complex 
topologies. [5, 11] The latter has been extensively studied on 2D surfaces. [12-16] 

In recent work we have extended and applied the methods of RC and DEC to derive a simplicial 
discretization of Hamilton's Ricci fiow [17]. As we apply RC and discrete forms to dynamic, but 
not necessarily covariant, geometric fiows on PF manifolds, it is necessary to develop a deeper 
understanding of the nature of curvature and exterior calculus in discrete geometries. In this 
manuscript we show how the curvature in PF manifolds gives rise to seemingly distinct notions 
of curvature tensors: (1) curvature with a single sectional curvature or (2) isotropic curvatures 
similar to those of Einstein spaces. Sec. I A and Sec. IB will review canonical RC and the standard 
approaches to discrete differential forms. Then in Sec. II we discuss the geometric principles 
behind hybrid cells as local measures and show how these hybrid volumes are core elements of a 
volume-based DEC. In Sec. Ill we discuss the representation of curvature operators over the hybrid 
measures and transformations between them. 

A. Canonical Regge Calculus 

Suppose A4 is a d-dimensional, smooth manifold endowed with a simplicial complex T. A PF 
triangulation of A^, To, is a mapping from each d-simplex to a flat d-simplex in IR'^ such that the 
proper lengths of the edges in the 1-skeleton of T are preserved in To- The simplicial manifold 
formed by To is often called a Regge manifold or Regge skeleton. We now review the canonical 
approach to RC, as in [6, 18]. 

The interior geometry of any given simplex in To is given by Euclidean or Minkowski geometry 
and represents a common tangent space for each of the vertexes of the simplex. Any such simplex 
has an induced metric that is uniquely determined by the proper squared-edge lengths of the 
simplex. The metric as a function of the edge lengths g^u{^'^) gives the local, piecewise-constant 
approximation to the metric associated with M. For two simplexes sharing a common boundary, 
the joint domain is isomorphic to and thus is intrinsically flat. While the two simplexes as viewed 
from an observer off the manifold may appear curved, i.e. with non-zero extrinsic curvature, a loop 



FIG. 1. 2-dimcnsional projection of a mapping of the neighborhood of a hinge h to R'' with defect Sh- 



of parallel transport from simplex A to simplex B and back induces no change in orientation on a 
tangent vector in general position. Formally this is related to the requirement of the existence of 
a metric compatible connection in the PF manifold. 

While there exists a flat connection across any (d— l)-boundary, curvature naturally arises when 
generating a map from a complete set of d-simplexes sharing a common (d — 2)-simplex to R*^. 
This is the first indication of curvature in the PF manifold. To effectively handle the discontinuity 
in mapping the neighborhood of a (d — 2)-simplex, or codimension-2 hinge h, we can map the 
neighborhood of /i to a subspace of M!^ and smoothly continue the mapping across the removed 
section of IR"^ resulting from 'breaking' a [d — l)-simplex into two. By breaking a (d — l)-simplex 
into two to make the mapping to R*^, a defect in the correspondence between the interior angles 
of the simplexes at h and 2n is evident (Figure 1). We take the deviation of the interior angles of 
each simplex on h from an exact embedding in flat space to be the defect angle £h associated to h 

£h = 2Tr-J20i. (1) 

This defect angle is the measure of curvature associated with parallel transport of a vector around 
a loop that encircles h. Indeed, if we take any vector with components in the plane h* orthogonal 
to h and transport it around the boundary of any area a'^^ (also with components in h*), then the 
vector will have rotated by an amount equal to e^, and the rotation occurs in the plane of h* . Any 
loop fj"'^ with components in h* will generate such rotations, independent of the area enclosed. This 
has implications for the sectional curvature associated with h. The sectional curvature associated 
with the loop a"^ is given by 

_ (Angle of Rot'n) _ Sh 
"'^ ~ (Area Enclosed) ~ {a'^f^l' ^ ' 

Since the angle of rotation is independent of the area enclosed, we can take an infinitesimal area of 
rotation encircling h whose limit, |cj"^| — t- 0, impling a singularity in the sectional curvature. The 
singularity is known as a conic singularity associated to the hinge h. 



We can provide a representation for the sectional curvature for a loop of parallel transport, but 
can we also provide a clear expression for the Riemann curvature? From the continuum theory 
with infinitesimal rotations we can express the rotation of a vector after transport along the 
boundary of an infinitesimal area by 



(3) 
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However, the curvature in PF manifolds is characterized by finite rotations. Friedberg and Lee [6] 
showed that a correspondence between the above notions of defect angle and a Riemann tensor 
can be approximately given by 

Ri2 uih) ^ ehS{xi)6{x2) « -Ri2 2i{h) (4a) 



Ri i{h) « eh6{xi)6{x2) « R2 2{h) (4b) 
R{h) « 2eh8{xi)5{x2) (4c) 

yMR{h) = 2eh5{xi)5{x2)^ (4d) 

where Eq. (4d) is exact and the coordinates {xi, X2} lie in the plane of h* . In this approach, where 
one smooths out the discontinuities, one generates a limiting sequence of surfaces approximating 
the plane h* and takes the limit to the PF surface. This generates Dirac delta distributions such 
that the curvature is evaluated only on the hinge h and zero elsewhere. While the curvature tensors 
are only approximate in this sense, the integrand of the Einstein-Hilbert action. 



is exact. Thus we can take the standard action principle and have an exact expression for the 
lattice geometry (first locally). 



= (5) 

where the integration is over a domain containing a single h. The global expression is obtained by 
summing over all hinges, 

^//MfiA = ^E^'.-4- (6) 



B. Discrete Exterior Calculus 



We now move away from the explicit geometry of PF manifolds to the properties of differential 
forms on simplicial manifolds. The framework we will follow below has stemmed from the work of 
Whitney [19] and, later, Bossavit [20]. Recent use of differential forms in PF manifolds corresponds 
to surface parameterization [11] or finite-element methods [4, 10]. The use of differential forms in 
the simplicial lattice has been independently used in RC (see for example [8, 21-23]). In the 
description below, we will follow some of the notation and conventions of [10]. 

The algebraic structure of the simplicial lattice is such that we have a natural representation of 
a discrete chain complex, 

\ ,W A .(-^-i) A ... A A A 0, (7) 
and a discrete co-chain complex, 

A 4 4 • • • 4 s('^-i) 4 4 0. (8) 



5 



Moreover, given a geometric dual lattice to the simplicial skeleton one can construct the dual chain 
and co-chain complexes. 

Given a p-form oj on the smooth manifold A4 , we seek a representation of uj on the PF manifold 
To- The approach taken in DEC discretizations is to take the smooth image of To on A^, i.e. T, and 
evaluate uj on the p-skeleton of T. The simplicial approximation to uj is obtained by integrating 
over individual elements of the p-skeleton of T and associating that value with the image of the 
element in To; 

(^uj\s^P^y.= ujo{s^P^) = J^^^uj. (9) 

Here we use the notation that denotes an element of the p-skeleton of T. In general, we 
will use s^'^ to denote both elements of T and To- In the standard RC approach to discrete 
differential forms, the discrete forms may often be denoted by subscripts, e.g. uje = {^\^) will 
generally denote a one-form uj associated with an edge i in the simplicial lattice. Only when we 
are considering the initial discretization will s^^^ be in T. Once the discrete differential form is 
assigned, all further manipulations take place within To- 

In addition to simplicial forms in To, one can also discretize uj on the dual lattice. Ambiguity 
arises here due to the multitude of ways of assigning a dual lattice. In general, there exist several 
natural geometric dual lattices, e.g. barycentric, circumcentric or incentric dual lattices, but there 
are also numerous non-intuitive dualities that can be constructed for an arbitrary simplicial lattice 
[24]. For our purposes and for clarity in later sections, we take the circumcentric dual as our 
prescribed dual lattice. The circumcentric dual of a /c-simplex is given by [25] 

^s^^) := sgn[Cfc,Cfc+i,...,Crf], (10) 

where Cp is the circumcenter of a and sgn is used to ensure an orientation consistent with 
the encompassing s^'^\ The circumcentric dual has several nice properties that make it a natural 
choice: (1) a simplicial element and its circumcentric dual are (locally) orthogonal to one another, 

(2) the p-dimensional dual to a p-element is equidistant from each vertex on the p-element, and 

(3) in special cases the circumcentric dual lattice corresponds to the Voronoi lattice generated by 
the 0-skeleton of the simplicial lattice. 

Once one has a dual lattice, discrete differential dual forms are obtained analogously to the 
simplicial forms. For a p-form ui and a p-element of the dual lattice, ui associated to C7(P) is 
given by 

^a;|>)) := w(>)) = J ^ ^u. (11) 

This provides a way to directly discretize continuous forms on the dual lattice. Our next step is 
to show how to transform discrete forms on one lattice to forms on the other. 

The discrete Hodge (or -k) dual is an isomorphism between the simplicial A;-forms and the dual 
(d — A;)-forms. Given that the definitions of the dual/simplicial forms are integrated quantities over 
their respective lattice elements, we can assign an average scalar density to the lattice element by 
dividing by the volume of the lattice element. The Hodge dual isomorphism then says that these 
average scalar densities are equal for two dual elements; 

' W^)\ = -r^(*a;|..('=)\, (12) 



where A^(k) = \s^^'^ \ denotes the integrated measure, i.e. norm, ofs(^-). If we take the dual lattice 
to be the circumcentric dual, then the relation in Eq. (12) says that integration over a volume 
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spanned by s^'^^and •ks^''^ is preserved under the Hodge dual, i.e. the discrete manifestation of the 
self-adjointness of -k in the L^-inner product. 

It is useful to note here a property of the Hodge dual in the lattice that will be useful later 
on. As in the continuum, if one scales a differential /c-form w by a scalar a, then the dual * (aw) 
has the same orientation as -ku. In the simplicial lattice the Hodge dual acting on a subspace of 
a lattice element s^'^^ returns back the lattice element -ks^^^ possibly with a scalar coefficient. In 
particular, if one takes the dual of the portion of a triangle t closest to an edge i, then one simply 
obtains the dual element to t itself. The dual lattice element t* orthogonal to t lies on the face of 
the polytope i* and so entirely lies in the convex hull of i and I*. Moreover it is orthogonal to t 
and therefore trivially orthogonal to the subspace of t closest to i. 

Given a general dual lattice to the simplicial lattice and the Hodge dual, we can construct the 
complex 

• • • 

containing the chain and co-chain complexes in the dual and simplicial lattices. In the above 
complex, the operators d and 5 define the co-boundary and boundary operators, respectively. 

The next ingredient in the calculus of discrete differential forms is the exterior derivative d 
which maps fc-forms to {k + l)-forms, 

s(fe)gs(fe+l) 

Here we have used the lattice boundary operator (via Stokes' Theorem) such that the action on a 
A:-element of the simplicial /c-skeleton returns its {k — l)-boundary. The discrete exterior derivative 
allows us to express dco in terms of the valuations of uj on the boundary of a given s^^^^\ By 
replacing the s^^\s^^^^^) with )) we can translate the discrete exterior derivative from 

simplicial forms to dual forms. 

Similarly the exterior coderivative, 6 = -kd-k, on a (/c — l)-form can be derived by using the 
adjoint relationship (in the continuum) between d and 6, 

= = ^ . (15) 

Again, the same relationship holds on the dual forms 

The co-derivative on discrete forms requires only our notion of the Hodge dual and the exterior 
derivative and so the above expressions come as a direct result of the adjoint relationships of these 
two operators in (-I-). We can, in fact, rebuild much of the algebra and calculus on exterior forms 
with these basic building blocks. The wedge product can also be reconstructed on the discrete forms 
as shown by Desbrun et al. [25]. We do not reproduce the wedge product here, but only mention 



(13) 
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that this wedge product will, in general, be non-associative in the discrete scale. Associativity, 
however, is recovered in the continuum limit. One can construct an associative wedge product, but 
such a construction will fail to preserve the anti-commutativity property. It is generally a feature 
of discrete physics with finite angles that certain continuum symmetries fail at the discrete level 
only to be regained in the infinitesimal edge-length limit. 

II. THE LOCAL STRUCTURE OF DISCRETE FORMS 

We have described thus far the previous approaches to analysis of curvature and differential 
forms on a To- What was evident when we examined the discrete forms and the relationship 
between forms on the dual lattices was that the integrated measure of a lattice differential form 
was preserved under the Hodge dual. Given this preservation, we can interpret the discrete forms 
as measures over a volume spanned by a lattice element and its dual, i.e. the convex hull of 
the vertexes of a lattice element and the vertexes of the dual element. This suggests a volume- 
based approach to discrete forms. Since many of the operations applied to discrete forms rely on 
transforming one discrete form on a fc-skeleton to a discrete form on a p-skeleton (where we may 
have p ^ k), we will be required to transform objects from one domain to a non-coinciding (but 
possibly overlapping) domain of integration. We now turn to the local properties of discrete forms 
and their inherent domains of support. This will allow us to reconstruct transformations between 
lattice elements in terms of hybrid simplicial-dual volumes, henceforth called hybrid cells or hybrid 
volumes. 



A. A Menagerie of Hybrid Cells 

A hybrid volume is a domain that is at once the local measure of a lattice element s*-'^^ and the 
orthogonal subspace dual to s^^\ For standard tensor analysis it is sometimes convenient to take 
the d-simplexes as the local domains over which a function or tensor is piecewise evaluated. This 
is due to the ability to define an unambiguous tangent space to any d-simplex in To- However, an 
arbitrary fc-form will, in general, be contained in multiple d-simplexes and thus requires junction 
conditions to hold across the multiple coordinate charts assigned to the d-simplexes sharing a given 
fe-simplex. We will show how hybrid cells can be viewed as natural domains for differential forms 
in To such that there are local orthogonal frames containing the carrier of the geometric content of 
the discrete form. The lattice element hybrid cells are atomized via local domains that contain the 
minimal, non-trivial amount of information about the lattice fe-skeleton. We build these hybrid 
cells from local constructions and discuss how they are representative of the measure of the discrete 
forms. 

We first construct domains that are shared by a set of fe-simplexes, one A;-simplex for each 
k = 0, . . . ,d. These shared domains become the "atoms" of our geometry in the sense that these 
are the simplest meaningful d-volumes in a PF manifold. Of course one could always subdivide 
these irreducible domains in some arbitrary way, even to go so far as to define infinitesimal domains. 
However, doing so yields no further discrete information about the lattice or the differential forms 
on the lattice. 

Definition 1. An irreducible hybrid cell, V^{o)^{i)...^{d) , in a PF d-manifold To is the d-simplex 

sgno...d([Co,Ci,...,Crf]) = -J-ei,...,,e^i A • • • A 
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FIG. 2. An irreducible domain in three dimensions. The domain (red, shaded) is that domain on which a 

0- form on A, a 1-form on [AB] ~ £, a. 2-form on [ASC] = t, and a 3-form on [ABCD] = T are mutually 
defined. The black edges bounding the irreducible domain are edges that lie in either the simplicial or dual 

1- skeleton of To- The points Ci, Ct, and Ct label the circumcenters of £, t, and T respectively. Moreover, 
there is a natural orthogonal basis in this domain given by the edges CyCi, CiCt, CtCr- This irreducible 
domain is foundation of the DEC as applied to the geometry of PF manifolds. 



given that s^*^^ G sC^^+i) every k < d, ei^i^-.-i^ is the orientation of s^'^\ and the vectors Cj are 
the vectors emanating from Cq. When a circumcenter lies outside To, we only take the domain of 
Vg(o)s(i)...s(d) that lies in To- 

Here a fc-simplex, or convex hull of /c + 1 points, is denoted [ao,ai,-- - , afc] and Ck is the 
circumcenter of the simplicial element s^^\ The factor sgUQ...^ ensures the volume is consistent 
with the induced orientation from s^'^\ i.e. if the orientation of [Co,Ci, . . . ,Cd] is opposite to that 
induced in s^'^^ then the orientation is flipped. This is matched by the totally-antisymmetric tensor 
^ii«2 - id whose orientation is induced by s^'^\ In cases where the circumcenter of an element lies 
outside the element, the volume gives negative contribution to sums over the irreducible hybrid cell. 
The 2-form still contains a vector with negative 1-dimensional orientation even after being made 
compatible with the containing volume's orientation. In general, the circumcenter lying outside 
the simplicial element leads to an over-counting of volume. This careful accounting of orientation 
ensures the conservation of total volume, e.g. over the simplicial element s^'^\ 

It is clear from the definition that the 0-skeleton of an irreducible hybrid domain consists of the 
circumcenters of each of the /s-simplexes that share the domain. In the 0-skeleton only the vertexes 
Co = s^^^ and = cr^^^ are vertexes that are also members of the simplicial or dual skeletons. 
Similarly, the 1-skeleton of Vs(o)g{i)...g(d) consists of 6 vectors, two of which are subspaces of either 
the simplicial or dual 1-skeleton of To- Figure 2 shows an irreducible hybrid cell in three dimensions 
and highlights the members of the 1-skeleton that lie in either of the 1-skeletons of To- 



In each irreducible hybrid cell we can form an orthogonal (or orthonormal) basis from the set 
of vectors {nii \mi = CjCj+i}- Moreover, since each of these cells are subspaces of R"^, we have the 
volume element given hy V = ^ejj...j^m*i A- • - Am*'*. In this basis of differential forms, we carry the 
information about the discrete 1-forms £ cx mg and A oc rrid-i- All other edges in the 0-skeleton of 
this cell are virtual in the sense that these 1-forms only become basis elements of higher-dimensional 
/c-forms while not carrying direct information about the discretization of 1-forms on the lattice. 
Similarly, the irreducible hybrid cells contain discretization content for only one simplicial and one 
dual fc-form, for all k. For simplicial fc-forms, this is clear from the definition. For dual A:-forms it 
is evident from the restriction of the cell to exactly one {d — A;)-simplicial element. The rest of the 
(fc+i) A;-forms in the irreducible hybrid cell carry the information about the graded algebra in the 
simplex s^'^\ but only indirectly given the flat interior of the simplex. 
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The irreducible hybrid domains give the smallest domain of support for any given k-foim in the 
dual lattices of To- From these irreducible domains we wish to reconstitute local "natural volumes" 
associated to each lattice element of To- The irreducible hybrid cells do tile the PF manifold To 
but will not generally provide a disjoint cover of To- Only when the triangulation is well-centered 
(i.e. when the circumcenter of each simplicial element is contained in the simplicial element) does 
the set of irreducible hybrid cells form a disjoint cover. We can now state a result with regard to 
the measure of the hybrid cell. 

Theorem 1. The convex hull (interior to To) of a simplicial element and its dual, CH(s^^\-ks^''^)n 
To, defines the domain of support for a discrete form on s^^^ (-ks^''^ ) which is given by the set 
theoretic union of the irreducible hybrid cells (contained in Tq). The volume measure of discrete 
forms on s^'^^ (s^''^ ) is given by Vg(k) = t^yIs'-'^^I | ^s'-'^^I- 

\k) 

Proof. For each irreducible hybrid cell containing a given s^^\ y...sik)..., any /c-form u will generally 
have a non-zero evaluation over s^^^ and hence have a non-zero component in each irreducible cell 
containing s^^\ We can check that the set union of K..^{fc)... is convex by examining the convex 
sum of extremal points on the ^..^(fc)...- First we examine the sum over s^'^ for p < k, which gives 

[J [Co, - - - Ck-i,Ck, . . . , Crf] n To = [vo,vi, - - - , Vk, Ck+i, . . . , Crf] n To, 



5(0) o(fc-l) 



s(fc) 



which is clearly convex. Consider first the convex sum of Cp and C'p for two s^^^s {p > k) for 
which the irreducible hybrid cell contain s^'^^ and s^^^ is non-zero. We can examine the two ir- 
reducible hybrid cells who only differ by s^p\ There exists a boundary between the two cells 
[vo, . . .,Vk,Ck+i, . . . ,Cp, . . . ,Cd] and [vq, . . . ,Vk,Ck+i, ...,Cp,.. . If p = d then the boundary is 
a subspace of the common s^'^~^\ In this case, the convex sum of the two Cd and C'^ is the dual edge 
A, which is a straight-line entirely contained on the combined domain (when the two d-simplexes 
are mapped onto R"^) . If p < d then we know from the circumcentric duality that the convex sum 



CpC' = pCp + (1 - p)C' p G [0, 1] 



forms angles ZCp^iCp,CpCp < ^ and ZCpCp,Cp+iCp < ^ and hence lies within both the irreducible 
hybrid cells. We therefore conclude that the set union of these irreducible cells for a given lattice 
element is CR{s^''\i^s^''^). 

However, this entire domain may not contribute to the final oriented sum over the V„g(k)..s. 
Two irreducible hybrid cells covering the same domain with opposite orientations will give zero 
contribution from the overlapping volume. When summing over the s^^-* for p < k, we simply get 
back a simplex 

^[Co, . . . ,Ck, ■ ■ ■ ,Cd\ = [ vo, ■ ■ ■ , Vk , Ck+i, ■ ■ ■ ,Cd- 

p<k 



s(fc) 



Meanwhile summing over the s^^ for p > k returns 

^[Co, ...,Ck,...,Cd] = [Co, ...,Ck, {ct(^-^^) }]. 

p>k 

Using these two results, the full sum gives a bipyramid with base a^^~^^ and /c-dimensional altitude 
given by s^''^ whose volume is Vk = -nx[s^''^[ | ★s^'^^j. □ 

\k) 



10 



The hybrid domains associated to lattice elements of To define the local measure for the discrete 
forms on the lattice. This domain defined by the measure may not actually encompass the lattice 
element, especially when the lattice is non-Pittway (when the element and its dual have empty 
intersection). However, the domain of support is generally more expansive and necessarily contains 
the lattice element. It is particularly insightful to notice that for a simplicial element and its dual 
the domains of support and local measures for a discrete form lo and its dual -koj coincide. These 
lattice-element hybrid measures and the convex domain of support then become fundamental to 
the discretization and algebra of discrete forms. Henceforth, when we think of the hybrid cell, we 
will take this to mean the local lattice measure and not the domain of compact support, since it 
is only the former that is necessary for explicit computations. Examples of the menagerie of the 
local measure hybrid cells in 3 dimensions is shown in Figure 3. 





Simplicial 


Dual 


Hybrid Cell 
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FIG. 3. In three dimensions there are four distinct classes of hybrid cells in To- We show in this figure 
representative depictions of these four classes. The labeling ((jj^^.^) labels the dimension of a A: simplicial 

elements'^'^) and its dual. The hybrid cell is heuristically constructed by connecting the vertexes of s'^^^ with 
the vertexes of the dual cell ★s^'^^. In the cases of (°) and (jj), the hybrid cell is just equal to s'^'^^ or a^'^\ 
respectively. 



We have built up the lattice element hybrid cells from irreducible domains. We now take a step 
in the reverse direction to examine hybrid cells common to multiple lattice elements, not just a 
given lattice element. Such constructions are very useful when we are required to relate discrete 
forms on a /c-skeleton to discrete forms from a p-skeleton {p 7^ k). For example, the exterior 
derivative duj is a map from fe-forms to {k + l)-forms and requires a relationship to be formed 
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between hybrid cells that overlap but do not coincide. The integration over the hybrid domain to 
gik+i) picks up contributions from each of the s^^^ € s^''~^^\ but they do not contribute equally. 
Rather the "democratic" allotment of domains by the hybrid cells for each s^'^) defines a domain 
of support for each of the s^'^^. The integration over the hybrid domain for the 5^^+^) then acts as 
a restriction on the integration for the s^*^^ and we find a domain common to s^'^^ and 5^^+^) for 
each term in the discrete dev. 

Corollary 1. The hybrid cell common to two simplicial elements s and s' (with dim{s) < dim{s')) 
is given by sum of oriented volumes 

+ X] K(o)...s-sW ^Km...s'-sW- (18a) 

s(fc)|^ s("i)|^, 

When the simplicial lattice is well-centered (i.e. c{s^'^^) € 5^*^^ Vs^'^^J, then the hybrid cell common 
to a set of simplicial elements S = Si^\ S2^\ . . . , Sn"^ is given by the sum; 

^l,...,n = X] K(o)---s{n)...s(in)...s(d)- (18b) 

Proof. The hybrid domain common to any two simplicial elements is obtained by the intersection of 
the hybrid cells for each element s and s'. This reduces the problem to pairwise set intersections over 
irreducible hybrid domains. For an irreducible hybrid cell containing both s and s' , the intersection 
returns the full irreducible cell. Any irreducible hybrid cell for s that is not also a hybrid cell for 
s' will yield the subspace common to both irreducible cells, but with an orientation opposite of 
the simplicial complex. These contribute with negative volume and remove subspaces not in any 
shared s^'^^ 3 s, s' . Hence, the hybrid cell common to s and s' is reduced to a sum over irreducible 
hybrid cells plus negative volume terms obtained by intersections of non-shared irreducible cells. 
The extension to n simplicial elements is straightforward, though quickly becomes cumbersome. In 
the case of well-centered simplicial complexes, all irreducible hybrid cells are disjoint and trivially 
factorize any given s^'^\ Hence, the second summation Eq. (18a) vanishes. □ 

The above definitions have only focused on hybrid cells for simplicial lattice elements. As we 
noted before for the hybrid cell for a given lattice element, the hybrid domain for s^^^ is coincident 
with the hybrid domain for ^s^'^^ as a result of the self-adjointness of the Hodge dual in the L^-inner 
product. Generalization from simplicial elements for the reduced hybrid cells or domains common 
to multiple lattice elements is, therefore, trivial by simply taking the dual of an element in the dual 
lattice. 

B. Solder Forms and Moment Arms 

We consider a frame bundle on our manifold. In each simplex we have a fibre that is a copy of 
the tangent space on the base manifold, which we call the space of values. The tangent space of 
the base manifold is the horizontal section while the fibre or space of values is the vertical section. 
On To we assign a tangent space to any d-simplex or any pair of neighboring d-simplexes. For each 
of these tangent spaces, we have a copy in the space of values. 

The PF manifold has a hard-wired simplicial skeleton and from that we construct a dual skele- 
ton, entirely determined by its rigid predecessor lattice. We can think of these two lattices as 
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complements of one another much the way we consider the space of vectors and one-forms as com- 
plements. Moreover, since the dual lattice is decomposed into subspaces of d-simplexes, we can 
define the dual lattice entirely in terms of the vectors in the tangent spaces. Thus, determining 
transformations between the simplicial and dual lattices is tantamount to determining an appropri- 
ate solder form in the d-simplexes, or more appropriately in the irreducible hybrid cells. A solder 
form, the unit vector-valued one-form, is the identity map from the tangent space~with one-form 
basis e'^-to the space of values—with vector basis Ca] 

dV = Cae". (19) 

When acting on a vector dV transforms a vector in the tangent space to a vector in the space of 
values. The summation ensures that we retain the basis expansion but in the new vector space. 
We can now construct a representation of the solder form in To that allows for transformations 
between the simplicial and dual lattices. 

In the discrete manifold, the solder form is an object that carries information about the rela- 
tionship between the space of values (vertical section) and the tangent space (horizontal section). 
In conjunction with the Hodge dual, the solder form acts as a transform between fe-forms, w, and 
p-forms {p < d — k) orthogonal to co. Through wedge products we can form a {d — p)-dimensional 
space orthogonal to a desired p-form that contains our /c-form. The Hodge dual and a summa- 
tion over discrete /c-forms that are able to construct such a space then provide a transformation 
between a set of /c-forms and the desired p-form. Another route is to take the Hodge dual of 
fc-forms orthogonal to a desired p-form and contract over the directions orthogonal to the p-form. 
Together these paths form two routes towards a notion of the double-dual of a vector-space valued 
differential (discrete) form. 

We start the construction in an irreducible hybrid cell. In this irreducible cell we have exactly 
one i and A that are representative edges of the simplicial and dual 1-skeletons, respectively. As 
discussed above, this domain has volume form given by 

Vs(o) ,e,- ,x,sW = ^(^31) ^ ^^^^ ^ ^) ' (20) 

where the M£x = (^z^^^i A • • • A 171^-2 is the {d — 2)-dimensional subspace orthogonal to both 
i and A, which we call the moment arm from £ to X. Both i and A are understood to consist of 
only the segment of i and A within a given irreducible hybrid cell. Taking the set {£, A, {mi}f~^} 
as an orthonormal basis (where {rui} are the one-forms that span Mix), we have the relationship 



between the lattice forms and their induced vectors; 

£{£) = 1, ^(A) = 0, £{mi) = (21a) 

A(A) = 1, A(^) = 0, X{mi) = (21b) 

miim) = 1, mii^) = 0, mi(A) = 0, (21c) 



where Mix represents any 1-form or vector in the subspace defined by Mix- We can then identify 
individual maps from each basis form to an identified basis element in the vertical section (space 
of values). In general, we make no real distinction between the lattice element as a scalar- valued 
differential form or as a vector-space valued differential form. In the context of this manuscript, 
the lattice elements always take the meaning of a vector-valued differential form, or as a map from 
the tangent space to the space of values. The solder form is thus 



"l(d-2) 

dPo = ^ + A(sW)+ mi(s('^)). 

mi 



(22) 
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The combination of the solder form and the Hodge dual provide us with practical tools for 
transforming between the two lattices. We first examine the use of the solder form in such a 
transformation that is well-known in general relativity. In E. Cartan's [2(3] approach to general 
relativity, the Einstein tensor need not be defined with relation to the Ricci curvature and Ricci 
scalar, but as the dual to the moment of rotation; 



G = * (dp A 7e) 



ieV^'^^/^^r'^rfSf, (23) 



where we have specifically worked in dimension 4. Here the solder form plays the crucial role of 
defining a subspace orthogonal to both the 3-volume for the moment of rotation and imposing the 
trace on the 3-form that results from the wedge product. In higher dimensions, we keep appending 
onto the Riemann curvature {d — 3) solder forms prior to taking the Hodge dual to obtain the 
moment of rotation {d — l)-form. In this light, the solder forms take on the interpretation as 
moment arms between the original basis elements and elements in the orthogonal subspace. This 
map provides a way to identify mappings between the dual lattices not provided by the Hodge dual 
alone. 

The solder form is an instrumental tool in the transformation of vector-space valued differential 
/c-forms to differential p-forms. Of particular interest to us now is the transformation of 1-forms 
uj G A(^) to 1-forms in the dual space uj' G A*(^\ which is the raising (or lowering) operation on 
differential forms. We take a lattice 1-form £ and construct a map to a dual one- form A. This is 
equivalent to mapping the 1-form £ to the dual space of 1-forms and asking for the components 
along A. 

To make such a map, we first make use of a discrete analog of a continuum property relating 
the Hodge dual, the wedge product A and the inner derivative i. 

Theorem 2. Let To he a PF simplicial manifold and choose a simplicial k-element s^'^^ G To, « 
one-form that extends from Ck to Ck+i of a s^^^^^ 3 s^^\ and a circumcentric Hodge dual 
operator, -k, on To- Then there exists a commutative diagram 



(24) 



using the wedge product with mk and the inner derivative (contraction) , Lm^, over the vector, m^, 
corresponding to ruk- 

Proof. We first prove this for the case k = \. Using the definition of the circumcentric Hodge dual 
and its property that any subspace of a p-simplex, s^^\ gives back -ks^^^ yields 

*(miA£)= ^ sgn[C2,C3,...,Crf] = ^ — m2 A • • • A mrf_i, (25a) 

s(p)B£,t: p>2 irii-. i>3 
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where the summations are over the s^^^ (for k > 3) containing i. Meanwhile, taking the Hodge 
dual first then the inner derivative results in 



f"mi (*^ 




nil A • • • A nid 



A rud 



(25b) 



ruj: j>2 



where we use the normalization property rrii (m^) = 5ij that picks out a given s^^^ from the 
summation. Multiplying by (d — 1) yields the desired result. 

The case for arbitrary k follows in a straightforward manner. □ 

This theorem shows that there exist two paths to map from an arbitrary /e-form on the simplicial 
lattice to a {d — k — l)-form on the dual lattice. These two paths stem from continuum descriptions 
of the double-dual transformations of tensors (as was used to construct the Einstein tensor). Since 
discrete forms are mapped as coefficients on the lattice elements, we need only know how the space 
of forms on the lattice elements get mapped to one another. The scalar coefficient gets carried 
through with no change. 

As a result of Theorem 2, we can show a subsequent commutative diagram for maps from A^^^ 
to A*(^\ We will drop combinatoric factors from the commutative diagrams for simplicity. 

Corollary 2. On To there exist maps from the simplicial 1- skeleton to the dual 1- skeleton such 
that 



Ad) 



> 



(26) 



p^*{d-i) ^\ ^*(i)^ 



Proof. The proof follows by induction and successive applications of Theorem 2 as in the diagram 
below: 



A(i) 



A(2) 



A(3) 



^*(d-3)^ 



Iteratively applying the above diagram, we finally construct the desired result. 



(27) 



□ 



The two paths in Eq. (26) define our two notions of the double-dual of a vector-space valued 
differential form. The Down-right path acts as the trace of the double dual to obtain a vector 
oriented in a given direction, while the Right-down path builds a (d — 1) subspace orthogonal 
to the desired vector before taking the dual to obtain the desired result. The Right-down path 
corresponds to the first line of Eq. (23) while the Down-right path corresponds to the last line of 
Eq. (23). 

Taking the inner-derivative of the hybrid volume form Vi\ with a desired A gives 



dixVix 



1 



(28) 



15 



which defines the orthogonal subspace to A. We find the moment arm by taking the inner-derivative 
again, 

d{d - l)LeLxVex = Mex- (29) 

The inner-derivatives are used to find the components of the moment arm that maps £ to A (or the 
reverse). Choosing another i or another A changes the moment arm. Using our above commutative 
diagrams, we then have 

A(i) A('^-i) 

(30) 

^*{d-i) ^ A*(i). 

Given that these moment arms play the same role as the solder forms in mapping A;-forms to p- 
forms, we will often refer to the moment arms as generalized solder forms. Just as the Einstein 
tensor can be expressed as 

G = * ( d V A-y Ad VATZ \ = Tr (*7^*) 

d— 3-times 

using the solder forms to map the 2-form to a 1-form, the moment arms A-Iix allow us to define a 
map from the simplicial 1-skeleton to the dual 1-skeleton on the dual lattice. In the more general 
case (Figure 4) of mapping a simplicial A:-form s to a dual p-form a {p < d — k), we define the 
moment arm Msa 

_ /d\/d-k\ d\ 

^-UJV P J'-'^^--^ k\{d-k-py.p& ''^^- 

This is the (p) -solder form between the simplicial and dual lattices. Of course, this solder form 
only makes sense when a and s have overlapping domains, i.e. when s € *cj or o" € Moreover, 
it is a natural consequence that the solder form between a lattice element and its dual is given 
by the point of intersection, and so we have a consistent framework to map between the dual and 
simplicial lattice elements. 

We can also construct solder forms between simplicial lattice elements. Given the duality 
between a and a simplicial element s^'^~^\ Eq. (31) also defines a moment arm between two 
simplicial elements. If we have two simplicial elements, s € {s^^^^} and s' G {s^P-*} then the moment 
arm or solder form between s and s' is given by 

Mss' = (^^ (^^ ~ isL^s'V^s,s'- (32) 



A similar result holds for solder forms between two dual lattice elements. This provides a 
geometric foundation for measures of discrete forms. The irreducible, reduced and standard hybrid 
cells form a topological framework on which we form local measures for lattice elements, while the 
moment arms/solder forms allow us freedom to transform from one lattice to another. This has 
yet to tell us anything about the structure of discrete differential forms in this context, and so we 
now shift our focus to the exterior calculus using the hybrid domains. 
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FIG. 4. Solder Forms are shown in three (left) and four (right) dimensions. When the dimensionahty of the 
space spanned by two elements s^*^-* and cr'^^ is (d — 1), a solder form between a s^*^) and a a^^^ is given by 
the vector from Cg(k) to C^g-i?) ■ more general cases, the generalized solder form between two elements is 
given by the (d — p — k) subspace orthogonal to both s^*^-* and ct'^-'. 

C. The Algebraic Structure of Forms on Hybrid Cells 

A common property of exterior calculus on discrete manifolds is the distributional nature of the 
discrete forms. In the canonical approach to DEC, we have discussed how the differential forms 
are obtained by integration of the continuous differential form over the simplicial element in T 
corresponding to a simplicial element in To- These discrete differential forms take on values when 
evaluated on the simplicial or dual skeletons. 

In applying a DEC formalism to curvature operators on the PF manifolds, we want to utilize 
the known properties of curvature in RC and use as our guide the Regge action principle. It is 
known from the canonical, continuum analysis of PF curvature that any loop of parallel transport 
in a plane orthogonal to a hinge will non-trivially transform tangent vectors carried around the 
loop. This occurs even when the loop of parallel transport does not intersect the orthogonal dual 
polygon to a hinge (independent of how one defines the dual lattice). As a result, the curvature 
remains non-zero as long as the loop of parallel transport contains a non-trivial projection onto a 
surface parallel to the dual polygon. Therefore, the discrete measure of the curvature in the hybrid 
domain of a hinge is given by 

[ RdVpropcr = ^ [ [ {mem ■ dh*) dh = 2ehAh, (33) 

JVh [2) Jh Jh* 

where the combinatoric normalization comes from the decomposition of the volume. We then 
notice that this is the total curvature across the hybrid cell. Yet it has little direct indication of 
the local tensorial content of the Riemann tensor. From a DEC perspective, we wish to find a 
projection of Riem into To that preserves this total curvature. 

Since the curvature in PF manifolds is entirely projected onto the polygonal dual h* to h we 
then require that this result be constant over h in order to recover Eq. (33). This is a direct 
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indication that the discrete form defines a constant field over h. If we take the Hodge dual on the 
discrete forms, then we also have, 



Riem • dh* 



dh = 2ehAh 



★Riem • dh 



The scalar coefficient is thus constant over the hybrid domain, and we now define an approach 
based on the standard DEC that explicitly assigns the discrete forms over the hybrid measures. 

In an explicit volume-based DEC, the base discretization is done as before by projection of the 
differential uj onto a simplex or dual polytope in T that corresponds to a simplex or dual polytope 
in %, 



[UJ\S 



{k)\ 



to (k) 



Our discrete forms then become the geometric objects ujg{k)S^^'^ where s^'^-^'s are the /c- forms in the 
hybrid cells for s^^\ The s^^'^ does not merely represent the lattice element, but the family of 
surfaces in Vg(k) parallel to s'^^\ each such surface mapped back to s^^'^ when viewed from outside 
Vg{k), i.e. when coarse-grained to smooth over the internal structure. As a comparison to the 
continuum, we have the relationship 



UJ, s 



(fe) 



ijj (k)dV, 



V 



using our previous definition of oo^^k). Then, using the projection of uj onto a lattice element, 
(wjs^'^)) = / ujg(k)S^'^'^ we have 



{k)\ ^ _±_ 
it) 



UJ s 



(34) 



We then take the coefficient uj^^k) as a scalar function defined over the dual polytope •ks^'^^ but with 
components only lying in the surfaces parallel to s^''^ The discrete measure of a /c-form uj is 



ui, s 



\kJ -^K^Cfe) 



(35) 



Further, this volume measure has the standard property that the Hodge dual preserves the coeffi- 
cient and so we have 

(★a;,*s(*^)) = (w,*(*s('=))) = (a;,s('^)) . (36) 

The definitions of the exterior derivative and co-derivative follow from Stokes' theorem applied to 
the L^-inner product; 



1 



(fc+i) •^K,(fc+i) 
1 1 



\k+l 



)k + 



(37) 
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where we have decomposed s('^+^) into s^^-^'s and the moment arms Mg{k)g(k+i)- Moreover, the 
integral is maintained over the measure corresponding to s^'^'^^^ and so the total volume integral is 
over the volume Vg(k)g(k+i). After integrating and combining appropriate terms, we then have 

\k+l) ^ ' s(fe)£s(fe+i) s(fc) 

Using similar properties of the co-derivative in the local inner-product, we obtain the measure of 
the discrete exterior co-derivative; 

(5a;, s^'^-^)) = (39) 

g(k)^g(k-i) 

Some key points and assumptions are useful to highlight before directly using this formalism. 
A core assumption for the DEC approach to vector-space valued differential forms is that we 
only explicitly discretize the base manifold. In this sense, we only discretize the tangent space 
components and examine the constraints this puts on components in the space of values. This 
is crucial for understanding the role of vector-space valued forms in the lattice. For a basis of 
the discrete forms, we expand a A:-chain in terms of the lattice elements, while the basis in the 
space of values is given by the unit p-forms with magnitude given by the coefficients defined in the 
discretization. Secondly, while discretization is done on a given lattice element, the /c-form field is 
extended as a constant field domain of support for the fe-element and measured on the hybrid cell. 
Under a standard assumption of the assignment of a field value to the lattice element, the discrete 
form would be distributionally valued on the lattice element. However, our measure requires the 
valuation to be constant across the orthogonal subspace to the lattice element. This is essential to 
properly recover local integral measures as averages over a given finite domain. 

We have provided a scheme based on standard DEC where the measures of /c-forms are given by 
integral d-measures instead of local /c-measures on /c-elements. This has the advantage of simplifying 
our evaluation of the discrete forms and providing a consistent framework for the understanding of 
the Hilbert action. However, we have only specified how one takes differential forms and discretizes 
them, but not how the inherent properties of the lattice affect the coefficients in the discretization. 
To examine the DEC approach to curvature in PF manifolds, we must examine properties of local 
curvature operators in addition to the integral measures provided by the Hilbert action. In the 
next section we analyze the standard curvature operators and the manifold geometry from the joint 
perspective of DEC and RC. 



III. CURVATURE FORMS IN PIECEWISE-FLAT MANIFOLDS 

Having laid out the foundations for analyzing curvature we can now provide measures of local 
components for the standard curvature operators. In particular, we examine the nature of curvature 
in the hybrid domain to codimension 2 hinges, simplicial edges, and simplicial vertexes. We have 
already discussed the role of the Hilbert action in measuring the curvature. We now specify how 
this can be obtained and used to derive local curvature operators. 

The Riemann tensor is defined as a (J) -tensor valued 2-form, or (by the raising operation in the 
space of values) as a bivector-valued 2-form; 

Riem = ^e^ A R^" ^^^e'" A = ^e^, A R^" „^ e" A e^. (40) 
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The Riemann tensor thus takes as an argument a bivector from the base manifold to return a 
rotation operator or rotation bivector in the space of values, i.e. for a loop of parallel transport 
enclosing a given area an arbitrary vector A = e^Af^ will be transformed as A — > A' = A + 5 A 
with 

In the discretization, the magnitude of 5 A will be determined by the geometry, while the basis 
bivectors in the space of values will remain a set of orthonormal bivectors. We can then set the 
basis of the tangent space as the integrated measures of the lattice elements. 

Similarly the Ricci curvature tensor is a vector-valued 1-form which is given by inserting an 
inverse solder form (a 1-form valued vector) into the bivector-valued curvature 2-form; 

Rc = trfp-iRiem = Riem(e^e/,) = e^R^'^^^e'^. (41) 

The inverse solder form induces a trace on the Riemann curvature tensor. Inserting the inverse 
solder form into the Ricci curvature yields the scalar curvature, 

R = i^-p-iL^p-iILiem = Riem(e'^e^, e^Cu) = 'Rc{e''e^) = R^" (42) 

In this section we will illustrate some properties of the representation of curvature via local 
operators in the PF manifold using the inherent discrete structure. Our general strategy is to 
formalize the notion of parallel transport and curvature interior to the hybrid cells at the scale 
lengths shorter than the local scale of the discretization. We then generate a representation of the 
Riemann curvature operator from the scale of the discretization, i.e. a zeroth order coarse-graining 
over the interior structure of the hybrid cells. This process is then used to reexamine the Ricci 
tensor from [271. 



A. The Riemann Tensor, Locally Einsteinian Structure, and the Conic Singularity 

Given the PF manifold To, the domain of support associated to a codimension 2 hinge, /i, is the 
convex hull (interior to To) of the hinge and a polygonal loop of parallel transport. Meanwhile the 
meaningful domain is the measure which is given by the oriented sum of irreducible hybrid cells. In 
an irreducible domain, there exists an orthogonal basis of vectors {mj} which are defined the same 
as before. Since the subspace of the hinge h is isomorphic to R'^"^, we can choose a complete set 
of the {rrii} for < i < d — 3 from any irreducible hybrid cell containing h. The last two vectors of 
any irreducible hybrid cell containing h then span h* . We take these two vectors as M^x and A. In 
the full measure of Vu-, we have a set of {Mh\, A} for each A G 5h* such that h* = ^ Z^AeSh* ^^^h\ 
(with a change in orientation as necessary to ensure consistent orientation across the domain). We 
instead focus on a basis of bivectors (and 2-forms) that leads to the representation of Riem as a 
(2) ^ (2) ™^trix. We then ask "How do we assign coefficients to the discretized Riemann tensor in 
this basis?" 

On the interior of the hybrid cell V/j the tangential components of the metric on the boundary 
between two neighboring tangent spaces are constant while there is generally a discontinuous change 
in the normal components of the metric across the boundary, when viewed in the basis of one of 
the d-simplexes [6, 28]. As we encircle a hinge with a loop of parallel transport, we notice that 
the components of the metric tangential to h are always constant across this domain (the surfaces 
parallel to h are always flat and remain parallel to h). Hence, the components of any vector in 
the space spanned by h will be unaffected by parallel transport around h. In our vector basis, 
the {rrii} for < i < d — 3 will be unaffected by parallel transport. However, the M^xs, which is 
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tangential to a given s^'^~^\ will generally have components normal to any other s^*^"^^ G Vh and 
so will experience a rotation when transported around a loop encircling h. The amount of this 
rotation is always given by the deficit, e/j. Moreover, this is true regardless of the size of the loop 
and depends only on whether the loop has a non-trivial projection into the plane orthogonal to 
h, i.e. h* . This indicates that the basis we have chosen is anomalous and that for a more robust 
analysis we must take care in our choice of basis. 

Theorem 3. The Riemann tensor, Riem, on codimension 2 hinges of To are rank 1 tensors with 
an eig en- decomposition 

Riem = r \'^\^ h*. 



V K 

Proof. We assign a basis {c*} such that for each we have 

a^.h* = 2g^^'^g'^Pal^^^hlp^Q (43a) 

where cr(^Q,) = \ [f^a + c^^] • Since each fi* has a non-zero projection onto h* , we have 

Riem(cJi) = h*£h. (43b) 

Therefore the Riemann tensor associates to each basis bivector Cj a rotation bivector oriented along 
h* with magnitude of rotation equal to e^- Here we have inserted an oriented area in Riem to 
obtain 

Riem(cj.j) = h*eh- (43c) 

This allows us to assign a matrix representation to the Riemann curvature tensor, 

( Eh {h*-a^) Eh {h*-a^) Eh{h*-a^) ■■■\ 
Eh {h* ■ a^) Eh {h* ■ a^) Eh{h* ■ a^) 
Eh {h* ■ ct'') Eh {h* ■ a^) Eh{h* ■ (7^) 

V : ; ; ■./ 



Riem 



(43d) 



where each row has identical elements and the sum of each column gives Ehh* . Since a Riemann 
tensor is a symmetric tensor across the basis of 2-forms, we require that the matrix representation 
adhere to the symmetry. The asymmetry that appears in Eq. (43d) is due to our asymmetric use 
of the conic singularity. To account for this, we further require that the inner product h ■ he 
normalized in both the space of values and tangent space. This is equivalent to requiring that any 
loop of parallel transport is treated as a 2-surface in the space with constant sectional curvature 
equal to the sectional curvature along h* . Hence any loop with non-trivial projection onto h* gets 
maximally projected onto the dual polygon h* as in Figure 5 and the Riemann tensor acts only on 
this projected loop. 



The Riemann tensor becomes a matrix with uniform entries, i.e. a rank 1 matrix with eigenvalue 
given by the sum over a row or column. Since there are (2) basis 2-forms, we have a Riemann 
tensor with a single component 



u d{d-l) ^ . d{d - 1) 
Riem = h — Eh h = — 



/ Eft • • • \ 







(43e) 
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FIG. 5. The conic singularity has the pecuhar property that the Riemann curvature tensor operates on any 
area of parahel transport (with non-zero projection onto h*) as if the area were projected identically onto 
h* . While we take the area of this projection to be that of ft,*, this is not a necessity of the PF manifold. 
All that is required is that each area only be acted upon by the sectional curvature of the plane h* . Taking 
the area to be equal to A/j. = is a normalization choice that is natural given the dual operation defined 
and the inherent orthogonality of the dual lattices, sec for example [17, 29]. 



If we then put the basis of the horizontal section as the full lattice elements h* and the orthogonal 
2-forms, then we have 

niem = h*^^^4^h*. (43f) 

□ 



In the eigenvalue analysis, we have assigned a projection operation on the hybrid cell that treats 
any non-trivial loop of parallel transport as a parallel transport around h*. This tells us that the 
Riemann tensor in this domain has the special property that there is exactly one non-zero sectional 
curvature in the plane h* , and a trivial flat subspace of bivectors orthogonal to h* . 

While this is now a local measure of the Riemann tensor as obtained from an arbitrary sampling 
of the space and evaluation of its eigenspace decomposition, it does not carry with it the flavor 
of the Riemann tensor as one would see from an evaluation on a single loop of parallel transport. 
This comes from the imposition of the conic singularity that all bivectors (with non-zero projection 
on h*) have assigned to them a equal sectional curvatures. Therefore, a basis in general position 
(all basis elements satisfying Eq. (43a)) acts as though the space were an Einstein space. The 
eigenvalue is the local measure over the hybrid cell and is the value through which we perform 
analysis in the DEC. 

From this eigen-decomposition of the Riemann tensor we can regain a local representation by 
normalizing the non-trivial eigenvalue by the combinatoric factor counting the basis elements of a 
basis in general position. In our case there are (2) 2-forms and the normalization factor is given 

by (2) = ^/{q}' This allows us to back-track from the eigenvalue ( as a measure of all surfaces 
in the domain) to the local evaluation of the Riemann tensor as a measure of the 2-surfaces in the 
domain of support; 

Riem = h* h*. 

Ah* 

This normalization process is a final step in the description of the discretized tensors on the PF 
manifold. It should be noted that any and all calculus is done on the local measures of the 
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differential forms and tensors, i.e. the tensors integrated over the local volume. The normalized, 
local tensors are convenient representations that are made possible by the locally simple structure 
of the hybrid cells. 

Further, we can assign a Ricci tensor and Ricci scalar in this hybrid domain. If we span h* by 
the local choice of {Mhx, A}, then we have 



which leads to Ricci tensor components 



d\ Eh 
2] Ah* ' 



and scalar curvature 



where the barred notation indicates the unnormalized representation of the eigenvalues. Normal- 
izing these curvature tensors given their differential forms character requires one to normalize by 
a factor of (2) for the Riemann tensor, for the Ricci tensor and (q) for the Ricci scalar. Hence 
the normalized curvature tensors are 

R^'^^.. =R\ = ^Il^ (45b) 

These are the oriented and normalized versions of the Riemann, Ricci and scalar curvatures within 
a hybrid domain Vh- We often only need the unoriented measure of the curvature forms and so 
summing over the orientations introduces a factor of 2 into Eqs. (44a), (44b), (45a) and (45b). 
We now have an accounting of the measures of curvature on the hybrid domain on a hinge that 
is analogous to that obtained by Friedberg and Lee [G]. On a qualitative scale, we have the same 
form of the Riemann, Ricci and scalar curvatures obtained in [(>]. The quantitative distinction 
comes from our use of the hybrid volume as a domain of support and our treatment of the Dirac 
distribution on h as spread out over the entire domain Vh, instead of distributionally valued only 
on h. 



B. The Ricci tensor and its double dual 

The Riemann tensor had a natural association to the dual polygons h* to the codimension 2 
hinges h given that its differential form properties are that of a bivector-valued 2-form. Hence 
one need only project two indices of the Riemann tensor onto the discretization. In evaluating the 
Riemann tensor in the simplicial discretization we also were able to express representations of the 
Ricci tensor and scalar curvature in this domain. We now shift attention to the PF representation 
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of the Ricci tensor. Since the Ricci tensor's natural representation is that of a vector-valued one- 
form, its direct discretization is on the 1-skeletons of the dual and simplicial lattices. In [27] we 
derived representations of the Ricci tensor in both the dual and simplicial 1-skeletons as weighted 
averages of the Riemann curvatures. We now try to elucidate the properties of these derivations 
and draw comparisons with an updated understanding of the curvature. 

The Ricci tensor is given by the bivector-valued 2-form curvature operator acting on the inverse 
solder form dT'~^ = e^e^, inducing a trace on the Riemann tensor. As a vector- valued 1-form, 
the Ricci tensor is directly discretized on 1-forms orthogonal to the hinges. Since Riem only has 
components in the planes h* , Rc only takes components along A or Mhx. The components along 
Mfix trivially give only one component from Riem and are contained in V/j. Calculating a Ricci 
tensor in the direction of M\h becomes a sum over directions orthogonal to Mxh, only one of 
which gives a non-zero contribution. Moreover, the M^x are virtual-being members of neither the 
simplicial or dual lattices- carry no inherent discrete differential forms. This is simply a statement 
about the non-independence of the curvature directed along M^x and the Riemann curvature 
associated to h. At the same time, in a small domain surrounding any given A, there are d distinct 
holonomies with independent curvature operators. Each of these distinct curvature tensors have 
Ricci curvature components oriented along A. Therefore, there exist non-trivial representation of 
Res, and distinct from the Riem of the hinges, on the A's of the dual lattice. In [27] we sought 
to ensure that the integrated measure of curvature associated h* and A was preserved over the 
domain common to these two elements. This is equivalent to the continuum requirement that 

tr (Rc, e'^)^ = tr (^Riem, e*" A e''^ ^ 

over some common domain fi. The measure of curvature one obtains from integration over the 
domain is given by the one non-trivial eigenvalue of the curvature in that domain, and hence the 
unnormalized measure. 

To trace the Riemann curvature in the domain of A is to sum over the polygonal loops h* 3 X 
given that the individual domains of overlap between A and each h* satisfies 

RxVxh' = Rh'Vh^x- (46) 
Doing so gives a scalar measure of the Ricci curvature on A 

Rx = (47) 

Given that this is an integrated measure, it samples all orientations of the loops of parallel transport 
and one naturally picks up the unnormalized, integrated form of the Riemann curvature tensor and 
an overall factor d{d — 1), 

R, = d{d-l)(^-^^^, (48) 
where we have used the volume-weighted average 

(A)a = • 

This is an association of a scalar quantity to a 1-form on the lattice. We can again normalize by 
the dimension of the space of 1-forms to obtain 

(49, 
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We can further associate to this measure a directionahty. As we have ah'eady mentioned, each 
of the curvature tensors contributing to R\ (Rx) is aheady oriented along A in each subdomain 
V\h since the measure of Rc on those domains is obtained by contraction of Riem with Mxh- 
Therefore, the measures R\ and Rx can be considered as coefficients on the one-form oriented 
along A. This assigns both directionality and magnitude to Rc on a given A. 

It is useful to note at this point that we have associated a scalar quantity to Rc on A and 
assigned to it a directionality. However, the scalar coefficient is dependent on the domain of 
integration. While we associate to A a component of Rc in the direction of A, this object is not 
of the same class as the Riem on the hinges. Whereas the curvature operators on the hinges are 
constant over the codimension 2 hinges and treated as constant over the domains Vh. (whenever a 
subdomain also encompasses the hinge itself), the Ricci curvatures are composed explictly in terms 
of components of curvature operators whose valuations only make sense within the subdomains 
Vxh. It is therefore notable that the scalar coefficient is only an appropriate measure whenever a 
domain of interest encompasses the entire V^- If a domain of interest only intersects a portion of 
Vx, then one must suitably restrict the measures in the definition of Rx (Rx)- This will come into 
play as we now seek representations of Rc on the simplicial lattice. 

The natural discretization of Rc is on the dual lattice; however, we have shown in [27] that a 
representation on the simplicial edges is also possible by requiring that domains of overlap between 
Vx and Vi give rise to the same measure of integrated curvature. We now want to show that is 
related to the discrete version of the 1-form double-dual of a vector-space valued 1-form. 

In each irreducible hybrid cell common to both A and i, we can form a basis from the two vectors 
i and A as well as the set of d — 2 vectors spanning the subspace M^x, {mi \ 1 < i < d — 2]. The 
transformation from the dual lattice to the simplicial lattice is done using the orthogonal subspace 
Mix applied through the commutative diagram from Eq. (26). If we take the discrete form along 
A, then we have {Rx) A as a vector- valued 1-form with A as the trivial map from A in the horizontal 
section to A in the vertical section. Taking the dual of this assigns Rx to A*. If a given (. is not 
contained in A* there is no contribution of Rx to Ri and likewise in the reverse. So we only consider 
when A € ^* and £ € A*. Moreover, Rx is constructed from the Riemann tensors evaluated on h*s 
and so we further expand the Rxs such that the dependence on fi/i* is explicit. If we then take 
the double dual and inner-derivative over the inverse solder forms on M^x-, we obtain 

RiVix = [iM,, i^Rx) , i) = {*Rx, i A Mix) = {Rx,*iiA Mex))y^^ = RxVex- (50) 

This object is a measure of the contribution from those curvature tensors with support in the defined 
domain. The above integrated curvatures are entirely connected to their domains of integration 
and the restriction of the domains induces a restriction of the integration on their definitions. 
We then view Eq. (50) as statement of dependence of the simplicial lattice Ricci curvature on the 
restriction of the dual skeleton Ricci curvature. As the latter depends on multiple curvature tensors 
in multiple domains, the restriction ensures that only those curvature operators with values in the 
specified domain contribute the final result. Summing over all A that are incident to i, i.e. all 
\ £ i*, gives the final measure of the unnormalized Ricci tensor on £; 

= Vi = vi = 'fe ^^'^ 



where 



and we have used the relation 



2^h3i ^hl 
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It is should be clear from Eq. (50) that the object being assigned to an i is not the component of 
Rc along i, but rather a measure of Rc in the orthogonal complement to i. Summing over all A's 
in i* provides a complete measure of the Rc in that orthogonal complement to i. 

After normalizing by the combinatorial factor (^) we get the normalized Ricci tensor on the 
simplicial edge £, 

R, = {d-l)^. (52) 

This demonstrates that we can formalize the transformation between the dual and simplicial lat- 
tices via the trace of the double dual. Moreover it provides a geometrically clear picture of the 
transformation via the moment arm or generalized solder form between ^ and A. 

The double-dual and its traces provide a direct path for the raising and lowering operators 
common in general relativity and differential geometry. These were applied in [27] as a method to 
obtain a dualized view of the Ricci tensor on the simplicial lattice. What is a particularly important 
lesson to be drawn from this is that the Ricci tensor associated with any I is not a measure of the 
Ricci tensor in the direction of I but an average of Ricci tensors in the {d— l)-subspace orthogonal 
to £. This was a crucial understanding in [17] since one must ensure that metric components in 
Hamilton's Ricci flow [30] change in proportion to those components of the Ricci tensor. If instead 
one were to simply take metric components along Hamilton's Ricci flow would not be recovered. 
Rather one would obtain an orthogonal flow to that of Hamilton's. 

We have done this explicitly for the Ricci curvature and similar results hold for the scalar 
curvature. In particular, taking the trace of the Ricci curvature, or the double trace of the Riemann, 
we assign a scalar curvature to a vertex, v, of the dual lattice; 

Since this is the scalar curvature, there is no distinction between the normalized and unnormalized 
curvatures and hence we drop the bars at the very beginning. In [•■Si] we showed that scalar 
curvature on a vertex v of the simplicial lattice is given by 

R, = d{d-l)j^. (54) 

We have thus described an intrinsic geometric derivation of the curvature operators and curvature 
scalar without explicit reference a limiting smooth sequence of surfaces. 



IV. DISCUSSION AND CONCLUSIONS 

We have shown in this manuscript a revised formulation of DEC that is based on the volume 
measures of differential forms on local domains of compact support. This formulation makes explicit 
the assignment of a discrete form to a family of surfaces in a volume local to a lattice element. 
The characterization of a discrete form to the family of surfaces provides a transparent view of the 
operations on discrete forms, such as the exterior (co)-derivative. This makes the DEC approach 
more directly amenable to use in RC. 

We built the volume-based DEC from the irreducible domains of the lattice-the monads of 
space-that form the most basic structures of the PF manifold. These irreducible domains are 
domains of supports for arbitrary fe-forms in any given tangent space. From these irreducible cells 
we identified generalized solder forms and local solder forms to allow for transformations between 
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the simplicial and dual lattices. Moreover, since solder forms provide a unit map from the tangent 
space (horizontal section) to the space of values (vertical section), a discretization of the solder 
form allows one to work within a framework that covers both scalar-valued and vector-space valued 
differential forms. 

We have also shown how curvature operators can be given explicit constructions as bivector- 
valued two- forms in the lattice by examining the conic singularities around the hinges /i of Tq. It was 
found that the PF Riemann curvature operators have an eigenspectrum with only one non-trivial 
eigenspace, that aligned along the plane orthogonal to the codimension 2 hinges. Noticing that 
the Riem operator takes a form analogous to a space with constant sectional curvature for bases 
in general position, we define local, normalized curvature operators that can then be compared 
to local continuum quantities. This provides flexibility to insist that not only should the average 
integrated curvature compare to the integrated curvature over a finite domain, but that the local 
curvature operators compare to sectional curvatures in the continuum. It was then shown how 
the Ricci curvature can be formulated on the dual lattice and how the simplicial representation is 
viewed not as the components of Rc along a simplicial edge H. but the components of the double 
dual of Rc. 

This general framework and the volume-based DEC is a purely discrete foundation for the anal- 
ysis of the geometry of PF manifolds. The continuum is only used at the level of the discretization 
and for the locally flat behavior of irreducible hybrid cells. We can therefore characterize this 
approach as a stratified view of discrete geometry with three distinct regimes: (1) the scale smaller 
than the local discretization where one admits ignorance of the internal structure and must assume 
some approximate behavior (e.g. flatness, constant curvature, etc), (2) the discrete scale where the 
geometric properties are based on the PF structure, and (3) a coarse-grained scale that regains the 
continuum behavior of the manifold. We have shown how to treat the discrete scale by inferring 
behavior from the connectivity of irreducible hybrid cells within a given domain. 

These results have been applied in [17] to a simplicial discretization of Hamilton's Ricci flow 
and in earlier stages in [32]. These results utilize the foundations laid in RC to open up DEC to a 
variety of geometric objects characterizable as vector-space valued differential forms. 
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